Coherent control in a decoherence-free subspace of a collective multi-level system 



M. Kiffner, J. Evers, and C. H. Keitel 
Max- Planck- Institut fiir Kemphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 

Decoherence-free subspaces (DFS) in systems of dipole-dipole interacting multi-level atoms are 
investigated theoretically. It is shown that the collective state space of two dipole-dipole interacting 
four-level atoms contains a four-dimensional DFS. We describe a method that allows to populate 
the antisymmetric states of the DFS by means of a laser field, without the need of a field gradient 
between the two atoms. We identify these antisymmetric states as long-lived entangled states. 
Further, we show that any single-qubit operation between two states of the DFS can be induced by 
means of a microwave field. Typical operation times of these qubit rotations can be significantly 
shorter than for a nuclear spin system. 
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INTRODUCTION 



The fields of quantum computation and quantum infor- 
mation processing have attracted a lot of attention due to 
their promising applications such as the speedup of clas- 
sical computations 0, 0, S| ■ Although the physical imple- 
mentation of basic quantum information processors has 
been achieved recently ^ , the realization of powerful and 
useable devices is still a challenging and as yet unresolved 
problem. A major difficulty arises from the interaction of 
a quantum system with its environment, which leads to 
decoherence [l, Q . One possible solution to this problem 
is provided by the concept of decoherence-free subspaces 
(DFS) 0, i, B E, [H [ll- Under certain conditions, a 
subspace of a physical system is decoupled from its en- 
vironment such that the dynamics within this subspace 
is purely unitary. Experimental realizations of DFS have 
been achieved with photons [13, [13, [3 [IB] and in nuclear 
spin systems [itI [isl . [l9| . A decoherence-free quantum 
memory for one qubit has been realized experimentally 
with two trapped ions [l^, [2l[ . 

The physical implementation of most quantum com- 
putation and quantum information schemes involves the 
generation of entanglement and the realization of quan- 
tum gates. It has been shown that dipole-dipole inter- 
acting systems are both a resource for entanglement and 
suitable candidates for the implementation of gate oper- 
ations between two qubits [13, [H, HI [H, HE [13, (Hi . 
The creation of entanglement in collective two-atom sys- 
tems is discussed in [12, [11] ■ Several schemes employ the 
dipole-dipole induced energy shifts of collective states to 
realize quantuni gates, for example, in systems of two 
atoms [M [H [ll, HI or quantum dots [li|. In or- 
der to ensure that the induced dynamics is fast as com- 
pared to decoherence processes, the dipole-dipole inter- 
action must be strong, and thus the distance between 
the particles must be small. On the other hand, it is well 
known that a system of particles which are closer to- 
gether than the relevant transition wavelength displays 
collective states which are immune against spontaneous 
emission [H, [1^ [30|, |3l|, [s^l ■ The space spanned by these 
subradiant states is an example for a DFS, and hence 



the question arises whether qubits and gate operations 
enabled by the coherent part of the dipole-dipole inter- 
action can be embedded into this DFS. In the simple 
model of a pair of interacting two- level systems, there 
exists only a single subradiant state. Larger DFS which 
are suitable for the storage and processing of quantum 
information can be found, e.g., in systems of many two- 
level systems [H, [s^l • 

Here, we pursue a different approach and consider a 
pair of dipole-dipole interacting multi-level atoms [see 
Fig. [1] . The level scheme of each of the atoms is modeled 
by a 5*0 <-+ Pi transition that can be found, e.g., in ^"Ca 
atoms. The excited state multiplet Pi consists of three 
Zeeman sublevels, and the ground state is a 5*0 singlet 
state. We consider arbitrary geometrical alignments of 
the atoms, i.e. the length and orientation of the vector 
R connecting the atoms can be freely adjusted. In this 
case, all Zeeman sublevels of the atomic multiplcts have 
to be taken into account [35| . Experimental studies of 
such systems have become feasible recently [36|, |37j, |38| . 

As our main results, we demonstrate that the state 
space of the two atoms contains a 4-dimensional DFS if 
the interatomic distance R approaches zero. A careful 
analysis of both the coherent and the incoherent dynam- 
ics reveals that the antisymmetric states of the DFS can 
be populated with a laser field, and that coherent dynam- 
ics can be induced within the DFS via an external static 
magnetic or a radio-frequency field. Finally, it is shown 
that the system can be prepared in long-lived entangled 
states. 

More specifically, all features of the collective two-atom 
system will be derived from the master equation for the 
two atoms which we discuss in Sec. |TT1 To set the stage, 
we prove the existence of the 4-dimensional DFS in the 
case of small interatomic distance R in Sec. IIIII 

Subsequent sections of this paper address the question 
whether this DFS can be employed to store and pro- 
cess quantum information. In a first step, we provide a 
detailed analysis of the coherent and incoherent system 
dynamics (Sec. IIV[) . The eigenstates and energies in the 
case where the Zeeman splitting 5 of the excited states 
vanishes are presented in Sec. lIV Al In Sec. lIVBl we cal- 
culate the decay rates of the collective two-atom states 
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which are formed by the coherent part of the dipole- 
dipole interaction. It is shown that spontaneous emission 
in the DFS is strongly suppressed if the distance between 
the atoms is small as compared to the wavelength of the 
So Pi transition. The full energy spectrum in the 
presence of a magnetic field is investigated in Sec. IIV CI 

The DFS is comprised of the collective ground state 
and three antisymmetric collective states. In Sec. |Vl we 
show that the antisymmetric states can be populated se- 
lectively by means of an external laser field. The proba- 
bility to find the system in a (pure) antisymmetric state is 
1/4 in steady state. In particular, the described method 
does not require a field gradient between the position of 
the two atoms. 

We then address coherent control within the DFS, 
and demonstrate that the coherent time evolution of two 
states in the DFS can be controlled via the Zeeman split- 
ting 6 of the excited states and therefore by means of 
an external magnetic field (Sec. IVI[) . Both static mag- 
netic fields and radio-frequency (RF) fields are consid- 
ered. The time evolution of the two states is visualized 
in the Bloch sphere picture. While a static magnetic field 
can only induce a limited dynamics, any single-qubit op- 
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FIG. 1: (a) The system under consideration is comprised of 
two atoms that are located at ri and r-i, respectively. The 
relative position R = r-i — r\ oi atom 2 with respect to atom 
1 is expressed in terms of spherical coordinates, (b) Internal 
level structure of atom £ {1, 2}. The ground state of each 
of the atoms is a So state, and the three excited levels are 
Zeeman sublevels of a Pi triplet. The states |2^) and |3^) 
correspond to the magnetic quantum numbers vrij = — 1, 
and 1, respectively. The frequency splitting of the upper levels 
is denoted by 5 = cjs — CJ2 = 0^2 — (j-"! , where hiOi is the energy 
of state 



eration can be performed by an RF field. 

In Sec. I VIII we determine the degree of entanglement of 
the symmetric and antisymmetric collective states which 
are formed by the coherent part of the dipole-dipole in- 
teraction. We employ the concurrence as a measure of 
entanglement and show that the symmetric and antisym- 
metric states are entangled. The degree of entanglement 
of the collective states is the same as in the case of two 
two-level atoms. But in contrast to a pair of two-level 
atoms, the symmetric and antisymmetric states of our 
system are not maximally entangled. 

A brief summary and discussion of our results is pro- 
vided in Sec. lVini 



II. EQUATION OF MOTION 

In the absence of laser fields, the system Hamiltonian 
is given by 

H = Ha+Hf+ i/vac , (1) 

where 

Ha = ?i^^a;.5|->|e, 

ks 

ffvac = ~d^'^ ■ E{ri) - df"^ ■ E{r2) . (2) 

In these equations. Ha describes the free evolution of the 
two identical atoms, huji is the energy of state and 
we choose Hcoi = 0. The raising and lowering operators 
on the \4fj} ^ transition of atom /i are {i G {1, 2, 3}) 

^1^ = M(4J and sl'^]^\A,){^,\. (3) 

Hp is the Hamiltonian of the unperturbed vacuum field 
and i?vac describes the interaction of the atom with the 
vacuum modes in dipole approximation. The electric 
field operator E is defined as 

ks 

where aks (ats) a-re the annihilation (creation) opera- 
tors that correspond to a field mode with wave vector 
k, polarization e^^ and frequency Wk, and V denotes the 
quantization volume. We determine the electric-dipole 
moment operator of atom fi via the Wigner-Eckart the- 
orem [3^ and arrive at 

d^"^ = ^ [d4^';> + H.c] , (5) 
1=1 
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where the dipole moments di = (*|d|4) are given by 
di = X'e(+) , d2= Ve^, 



(6) 



and V is the reduced dipole matrix element. Note that 
the dipole moments di do not depend on the index /i 
since we assumed that the atoms are identical. 

With the total Hamiltonian H in Eq. ^ we derive a 
master equation for the reduced atomic density opera- 
tor g. An involved calculation that employ s the Born- 
Markov approximation yields [H, [H, |4l| 

dtg^-^[HA,g]~^[Hn,g] + C^g. (7) 

The coherent evolution of the atomic states is deter- 
mined by Ha + Hq, where Ha is defined in Eq. ([2]). 
The Hamiltonian Hq arises from the vacuum-mediated 
dipole-dipole interaction between the two atoms and is 
given by 



Hn = -h 



3 

E 

i=l 



H.c.j 



0(1) 0(2) 



-h{n32 ( 



c-(2) oil) 
'-'3 + '-'2 - 



c.(2) 



H.c.j 
H.c.j 

H.c.j . (8) 



The coefficients ftij cause an energy shift of the collective 
atomic levels (see Sec. IIV[) and are defined as [H, HO, lU 



dj Re(x) d* 



(9) 



Here x is a tensor whose components Xki ^'^^ ^ ^ 
{1,2,3} are given by 



XkiiR) 



h.3 



dki \ - + ^ 7 

1] 77^ T]-^ 



RkRi fi 



i?2 



3i 



(10) 



R denotes the relative coordinates of atom 2 with respect 
to atom 1 (see Fig. [IJ, and 77 = koR. In the derivation of 
Eq. ([To]), the three transition frequencies wi, 012 and 073 
have been approximated by their mean value ujq = c/cq (c: 
speed of light). This is justified since the Zeeman split- 
ting S is small as compared to the resonance frequencies 
LUi. For i = j, the coupling constants in Eq. ([9]) account 
for the coherent interaction between a dipole of one of the 
atoms and the corresponding dipole of the other atom. 
Since the 3 dipoles of the system depicted in Fig. [T{b) are 
mutually orthogonal [see Eq. ([6])] , the terms f2y for i ^ j 
refiect the interaction between orthogonal dipoles of dif- 
ferent atoms. The physical origin of these cross-coupling 
terms has been explained in [4fl|]. 



The last term in Eq. ([7]) accounts for spontaneous emis- 
sion and reads 



2 3 



fi—li—l 

- E {r^^ + ^S^SI^ - 2S^gS^ + H.c.j 

i=l 

+ Far (81^18^^13 + gSi^S^^l 2S^^lgS'i'l 
+ F32 {Siisi% + gSiisi^l - 2Si%sii 

•}■ 



H.c. 



(11) 



The total decay rate of the exited state |i) of each of the 
atoms is given by 2^i, where 



7i 



1 2\d,\^L0^ 



7: 



(12) 



and we again employed the approximation uji ^ luq. 
The collective decay rates Fy result from the vacuum- 
mediated dipole-dipole coupling between the two atoms 
and are determined by 



di Im(x) d* 



(13) 



The parameters F^; arise from the interaction between a 
dipole of one of the atoms and the corresponding dipole 
of the other atom, and the cross-decay rates F.y for i ^ j 
originate from the interaction between orthogonal dipoles 
of different atoms [i^ . 

In order to evaluate the expressions for the various 
coupling terms f2y and the decay rates F^ in Eqs. ^ and 
(|13p . we express the relative position of the two atoms in 
spherical coordinates (see Fig. [1]), 

R — R (sin 9 cos 0, sin sin 0, cos 0) . (14) 

Together with Eqs. (fTU]) and ^ we obtain 



^^31 = 7^~3 [{^'^ ~ cos?/ — 3?/ sin 77] sin^ 9e , 



^21 = 

^22 

^32 



\(3ri^ - 1 + (r]2 - 3) cos 26*) cos?? 



— ?/ ( 1 + 3 cos 20) sin rj] 



-\/2 cot r^sie 



i4> 



nil - (2cot2 6l- 1)^2316^*' 
— ^21 , i^33 = ^^11 , 



(15) 
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and the collective decay rates are found to be 
3 

Tgi = 7— [{if - 3) sinr/ + Sr^cosTy] sin^ Oe^^"^ , 

Til = 3^ [{3rf - 1 + (772 - 3) cos 26*) sin 77 

+f] (1 + 3 cos 29) cos 77] , 

T21 = -\/2cot6ir3ie*'^, 

r22 =rn - (2cot2 0- I)r3ie2^'^, 

= -r2i , Taa = Til . (16) 

The coupling terms f2ii, Q31 and the collective decay 
rates Fn, T31 are shown in Fig. [2] as a function of the 
interatomic distance R. 

Finally, we consider the case where the two atoms are 
driven by an external laser field, 

El = [£a + SyBy] e'''^ '- e-'^^' + c.c. , (17) 

where S^, £y and e^,. By denote the field amplitudes 
and polarization vectors, respectively, is the laser fre- 
quency and c.c. stands for the complex conjugate. The 
wave vector — k^Ez of the laser field points in the 
positive z-direction. In the presence of the laser field and 
in a frame rotating with the laser frequency, the master 
equation (O becomes 

dtg = -^[Hl + Ha, g] - '-[Hn, g] + C^~g. (18) 

In this equation, is the transformed Hamiltonian of 
the free atomic evolution, 
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FIG. 2: (a) Plot of the vacuum-induced coupling terms Sin 
and according to Eq. psp . Ao is the mean transition 
wavelength. If the interatomic distance R approaches zero, 
the parameters iln and fiai diverge, (b) Plot of the collective 
decay rates Fn and Fai according to Eq. (|16p . Fn and Fai 
remain finite in the limit i? — > 0. The parameters in (a) and 
(b) are given hy 9 = 11/2 and = 0. 



HA = ~nj2jZ^^S['^S^r}. (19) 
i—1 /i— 1 

The detunings with the state \i) are labeled by A; = ujl — 
uii {i G {1, 2, 3}), and we have Ai = A2 + S, A3 — A2 — S. 
The Hamiltonian Hl describes the atom-laser interaction 
in the electric-dipole and rotating-wave approximation, 

2 

+ [-Q,ir,,) + tny{r^)] 5^^^ + H.c.} , (20) 

and the position-dependent Rabi frequencies are defined 
as 

n^{r) = V£J{V2h) exp[ifci -r] , 

ny{r) = V£y/{V2h) exp[ikL-r] . (21) 



III. DECOHERENCE-FREE SUBSPACE 

In this section we show that the system depicted in 
Fig. [1] exhibits a decoherence-free subspace. By defini- 
tion, a subspace V of a Hilbert space H is said to be 
decoherence-free if the time evolution inside V is purely 
unitary For the moment, we assume that the 

system initially is prepared in a pure or mixed state in 
the subspace V. The system state is then represented 
by a positive semi-definite Hermitian density operator 
gv e End(V) with Tr{gv) = 1- It follows that V is 
a decoherence-free subspace if two conditions are met. 
First, the time evolution of g^ can only be unitary if the 
decohering dynamics is zero, and therefore we must have 

C^gv = (22) 

for all density operators gv that represent a physical sys- 
tem over V. Second, the unitary time evolution governed 
by Ha + i?n must not couple states in V to any states 
outside of V. Consequently, V has to be invariant under 
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the action of Ha + Hci , 



e V 



{HA + Hn)\i') e V. 



(23) 



Note that since {H^ + Hq) is Hermitian, this condition 
also imphes that it cannot couple states outside of V to 
states in V. 

In a first step we seek a solution of Eq. ((22l) . To this end 



we denote the state space of the two atoms by Tisys and 
choose the 16 vectors = |ji)®|j2) {i, j G {1,2,3,4}) 
as a basis of Tisys- The density operator g can then be 
expanded in terms of the 256 operators 

\i,j){kj\, I, J, fc, {1,2,3,4}, (24) 

that constitute a basis in the space of all operators acting 
on Hsvs, 



4 4 

XI XI 

i,j=l k,l=l 



(25) 



It follows that g can be regarded as a vector with 256 
components gij^ki and the linear superoperator C-y is rep- 
resented by a 256 x 256 matrix. Equation ([22| can thus be 
transformed into a homogeneous system of linear equa- 
tions which can be solved by standard methods. 

For a finite distance of the two atoms, the only exact 
solution of Eq. is given by |4,4)(4,4|, i.e. only the 
state |4,4) where each of the atoms occupies its ground 
state is immune against spontaneous emission. A dif- 
ferent situation arises if the interatomic distance R ap- 
proaches zero. In this case, the collective decay rates 
obey the relations 

lim r^i = lim = lim = 

R~*o R^o R^o 

lim Fii = lim T22 — lim Ess = 7 . (26) 

fl^O R-^Q R^O 



In order to characterize the general solution of Eq. 
in the limit i? — > 0, we introduce the three antisymmetric 
states 



K,4) -|4,z)] , {1,2,3}, 



as well as the 4 dimensional subspace 

V = Span(|4,4), |ai), jaa), hs)) ■ 



(27) 



(28) 



The set of operators acting on V forms the 16 dimensional 
operator subspace End(V). We find that the solution of 
Eq. in the limit i? ^ is determined by 







6 e End(V) , 



(29) 



In particular, any positive semi-definite Hermitian oper- 
ator g\> S End(V) that represents a state over V does not 
decay by spontaneous emission provided that i? — > 0. 

We now turn to the case of imperfect initialization, i.e., 
the initial state is not entirely contained in the subspace 



V. Thenj_states outside of V spontaneously decay into 
the DFS [I4] . This strictly speaking disturbs the unitary 
time evolution inside the DFS, but does not mean that 
population leaks out of the DFS. Also, this perturbing 
decay into the DFS only occurs on a short timescale on 
the order of at the beginning of the time evolution. 

These results can be understood as follows. In the 
Dicke model t23(, i32,] of two nearby 2-level atoms, the an- 
tisymmetric collective state is radiatively stable if the in- 
teratomic distance approaches zero. In the system shown 
in Fig. [Tl each of the three allowed dipole transitions in 
one of the atoms and the corresponding transition in the 
other atom form a system that can be thought of as two 
2-level atoms. This picture is supported by the fact that 
the cross-decay rates originating from the interaction be- 
tween orthogonal dipoles of different atoms vanish as R 
approaches zero [see Eq. ([26])]. Consequently, the sup- 
pressed decay of one of the antisymmetric states \ai) is 
independent of the other states. 

In contrast to the cross-decay rates, the coherent 
dipole-dipole interaction between orthogonal dipoles of 
different atoms is not negligible as R goes to zero. It is 
thus important to verify condition (|23p that requires V to 
be invariant under the action of Ha + Hq . To show that 
Eq. holds, we calculate the matrix representation 
of Hfi in the subspace A spanned by the antisymmetric 
states {|ai), |a2), Ws)}, 



^^11 ^^21 ^"^31 

[HnU = ^ I ^21 ^22 

^31 ^32 ^33 



(30) 



Similarly, we introduce the symmetric states 

|s,) = -i=[|z,4) + |4,i)], ze {1,2,3}, (31) 

and the representation of Hq on the subspace S spanned 
by the states {|si), |s2), jsa)} is described by 



(32) 







0* 
"21 


fill 


[Hn]s = -h 1 


^21 


^22 


^32 




{ f^31 


^32 


f^33 



It is found that Hq can be written as 
3 

Hn = {ai\Hn\aj)\at){aj 

3 

+ V {s,\Hn\s,)\s,){s,\ 



(33) 



i.e., all matrix elements {ai\H^\sj) between a symmetric 
and an antisymmetric state vanish. This result implies 
that Hq couples the antisymmetric states among them- 
selves, but none of them is coupled to a state outside of 
A. Moreover, the ground state |4, 4) is not coupled to 
any other state by Hq. It follows that the subspace V is 
invariant under the action of 
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It remains to demonstrate that V is invariant under 
the action of the free Hamihonian Ha in Eq. With 
the help of the definitions of \ai) and \si) in Eqs. ((27|) 
and pip , it is easy to verify that Ha is diagonal within 
the subspaces A and S. In particular, Ha does not in- 
troduce a coupling between the states \ai) and \si), 



s,\HA\a,) = - [{iA\HAVA) - (4,i|iJA|4,i)] 



(1) c.(i) 



0. 



(34) 



Note that these matrix elements vanish since we assumed 
that the two atoms are identical, i.e. we suppose that the 
energy hwi of the internal state does not depend on 
the index ^ which labels the atoms. 

In conclusion, we have shown that the system of 
two nearby four-level atoms exhibits a four-dimensional 
decoherence-free subspace V C Ti-sys if the interatomic 
distance R approaches zero. However, in any real situa- 
tion the distance between the two atoms remains finite. 
In this case, condition Eq. (|22p holds approximately and 
spontaneous emission in V is suppressed as long as R is 
sufficiently small. In Sec. lIVBi we demonstrate that the 
decay rates of states in V are smaller than in the single- 
atom case provided that R < 0.43 x Aq. 



IV. SYSTEM DYNAMICS-EIGENVALUES AND 
DECAY RATES 

The aim of this section is to determine the energies and 
decay rates of the eigenstates of the system Hamiltonian 
Ha + Ha- In a first step (Sec. HV^ . we determine the 
eigenstates and eigenvalues of H^,- It will turn out that 
these eigenstates are also eigenstates of Ha, provided 
that the Zeeman splitting of the excited states vanishes 
{5 = 0). Section HV Bl discusses the spontaneous decay 
rates of the eigenstates of and Scc. lIVCl is concerned 
with the full diagonalization of Ha + Hn for 5^0. 



A. Diagonalization of Ha 

We find the eigenstates and eigenenergies of Hu by the 
diagonalization of the two 3x3 matrices [Ha]^ and [Hq\s 
which are defined in Eq. (|5D|) and Eq. ([5^ . respectively. 
The eigenstates of Hq, in the subspace A spanned by the 
antisymmetric states are given by 



where 



sin 0\a2) 
cos9\a2) 



C0s6l|V'a), 



-sin 6*1-0^ ) , 



V2 



[e'^|ai)±e-^*|a3): 



(35) 



(36) 




FIG. 3: Plot of the vacuum induced energy shifts Qf and 
n AT as a function of the interatomic distance 7? according to 
Eq. (|38|l . These shifts enter the expressions for the eigenvalues 
of Hst in Eqs. dST)) and (gT}. Note that ^Ip decreases with 1/R 
for large values of R, while Qn vanishes with \/R? . 



We denote the eigenvalue of the state |V'a) by and find 



A„ — A„ 



(37) 



where 



V ~ + ?7sin(77)] 



= 7^ [cos(77) -f 7/sin(77)] , 



(38) 



and rj = k^R. The parameters JIf and Qat are shown in 
Fig. [3] as a function of the interatomic distance R. 

The eigenstates of Hq, in the subspace S spanned by 
the symmetric states are found to be 



IV's) = sin 6*152) — cos 0|?/'s ) 
= cos 6*152) + sin 6*1 



where 



1 



\i^f)^-j=[e^^\s,)±e-^^\s,)] 

and the corresponding eigenvalues read 

Al A? = -;ir2p, Xl^-Kl 



N ■ 



(39) 



(40) 



(41) 



Next we discuss several features of the eigenstates and 
eigenenergies of Hq. First, note that two of the sym- 
metric (antisymmetric) states are degenerate. Second, 
we point out that the matrices [Hq\a and [-ffojs consist 
of the coupling terms f2y which depend on the inter- 
atomic distance R and the angles 9 and cf) [see Fig. [1] and 
Eq. (fT5|) ]. On the contrary, the eigenstates \tpD and \4>l) 
depend only on the angles 9 and </), but not on the in- 
teratomic distance R. Conversely, the eigenvalues of H^ 
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are only functions of the atomic separation R and do not 
depend on the angles 9 and 0. This remarkable result 
is consistent with a general theorem (s^l that has been 
derived for two dipole-dipole interacting atoms. The the- 
orem states that the dipole-dipole induced energy shifts 
between collective two-atom states depend on the length 
of the vector connecting the atoms, but not on its ori- 
entation, provided that the level scheme of each atom is 
modelled by complete sets of angular momentum multi- 
plets. Since we take all magnetic sublevels of the Sq <-> Pi 
transition into account, the theorem applies to the sys- 
tem shown in Fig. [TJ 

In Sec. llVCi we show that the eigenstates ItpD and 
Itpl) of Hn are also eigenstates of Ha, provided that the 
Zeeman splitting S of the excited states vanishes. This 
implies that the energy levels of the degenerate system 
{S — 0) do not depend on the angles 9 and (j), but only 
on the interatomic distance R. From a physical point of 
view, this result can be understood as follows. In the ab- 
sence of a magnetic field (6 — 0), there is no distinguished 
direction in space. Since the vacuum is isotropic in free 
space, one expects that the energy levels of the system 
are invariant under rotations of the separation vector R. 
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FIG. 4: Dependence of the parameters FJj and F* on the 
interatomic distance R according to Eq. (|43p . (a) In the limit 
7? — » 0, the F* tend to zero, and the antisymmetric states ) 
are subradiant. (b) The symmetric states Itpl) decay twice as 
fast as compared to two independent atoms if R approaches 
zero. 
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FIG. 5: The atoms are aligned in a plane spanned by the unit 
vectors and — (cos sin 0, 0). Within this plane, the 
relative position of the two atoms R — z e^ + l is described 
by the parameters z and I. The energies of the eigenstates of 
Ha + Hn depend only on z and /, but not on <j). 



B. Decay rates 

In order to find the decay rates that correspond to 
the Eigenstates j?/"*) and \tpl) of the Hamiltonian Hq, we 
project Eq. ITTI) onto these states and arrive at 

dtmgm = -2ri {i;l\g\^Pl)+Cl{t). (42) 

In these equations, 2 F^ and 2 F* denote the decay rates 
of the states and respectively. The time- 

dependent functions C^(i) and Cl{t) describe the in- 
crease of the populations (V'alfi'lV'o) {''PsIqI'^I) due 
to spontaneous emission from states \i,j) G {1, 2, 3}) 
where both atoms occupy an excited state. The explicit 
expressions for the coefficients F^ and F* as a function of 
the parameter rj = k^R are given by 

Tl = F2=7^[2r,3-3?7COs(r,)+3(l-772)sin(77)] , 

Tl = 7^ [r;3 + 37ycos(ry)-3sin(7])] , 

Tl = F2=7^[2773 + 377COs(77)-3(l-772)sin(7y)] , 

Tl = [rf - 377COs(r,) 3sin(77)] . (43) 

These functions do not depend on the angles 9 and (/), but 
only on the interatomic distance R. As for the dipole- 
dipole induced energy shifts of the states and \'4'\) 
(see Sec. IIV A[) . this result is in agreement with the the- 
orem derived in [35] . 

Figure ID^a) shows the parameters F^ as a function of 
R. The oscillations of F^ and F^ around 7 are damped 
with 1/R as R increases, and those of F^ decrease with 
1/R^. Note that the oscillations of the frequency shifts 
display similar features for R ^ Xq (see Sec. IIV Al) . 
It has been shown in Sec. IIIII that any state within the 
subspace A of antisymmetric states is completely stable 
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FIG. 6: (Color online) Plot of the energy shifts that determine the energy levels of the antisymmetric states according to 
Eq. (|45p . In (a)-(c), the parameters are shown in a plane spanned by and = (cos sin <jl), 0). The relative position 
i? = 2 + Z 60 of the atoms in this plane is parameterized by z and I (see also Fig. [5]). Since the AJ, do not depend on (j), the 
energy surfaces shown in (a)-(c) do not change if is rotated around the z-axis. While A^ and A^ tend to -co in the limit 
i? ^ 0, A^ tends to +oo. The frequency splitting of the excited states is 5 = 7. In (d), the A^ are shown as a function of the 
interatomic distance R, the parameters are 6 = -k 12 and 5 = 7. 



for i? ^ 0. Consequently, the decay rates 2rjj of the 
states IV'a) tend to zero as R approaches zero. It can be 
verified by numerical methods that V\ and Y\ are smaller 
than the parameter 7 provided that R < 0.44 x Aq, and 
does not exceed 7 if i? < 0.72 x Aq. For i? = 0.1 x Aq, 
the coefficients T\ are smaller than 0.1 x 7. Although R is 
larger than zero in an experiment, the states |^^) decay 
much slower as compared to two non-interacting atoms 
if R is sufficiently small. This shows that spontaneous 
emission can be strongly suppressed within the subspace 
A of the antisymmetric states, even for a realistic value 
of the interatomic distance R. 

The parameters Y\ are depicted in Fig. Sfb). In the 
limit R ^ 0, the coefficients T^ tend to 27. The sym- 
metric states within the subspace S display thus super- 
radiant features since they decay faster as compared to 



two independent atoms. 



C. Non-degenerate System 

Here we discuss the diagonalization of Ha + Hq in the 
most general case where the Zeeman splitting 5 of the 
excited states is different from zero. The matrix repre- 
sentation of this Hamiltonian with respect to the states 
{|^a), li'l), defined in Eq. m reads 

(ivo + np S cos 6 \ 
ScosO ujo + ^F —Ssm9 
—Ssm6 luq + Qn J 

(44) 
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In general, the eigenvalues of this matrix can be written 
in the form 

El - h{^o + Al) , 

El = ri{ujo + Al) , (45) 

where the frequency shifts depend only on the inter- 
atomic distance R and the azimuthal angle 9, but not 
on the angle 0. To illustrate this result, we consider 
a plane spanned by and = (cos </>, sin (/>, 0) , see 
Fig. [5] Within this plane, the vector R — z Ez + le^ 
is described by the parameters z and I, and Fig. [B{a)-(c) 
shows AJ,(Z, z) as a function of these variables. Since the 
do not depend on (j), the energy surfaces shown in 
Fig. [51[a)-(c) remain the same if is rotated around the 
z-axis. This result follows from the fact that the Hamil- 
tonian Ha in Eq. ([2]) is invariant under rotations around 
the z- axis [35|. 

In Sec. lVIi we will focus on the geometrical setup where 
the atoms are aligned in the x- y- plane (0 = 7r/2). In 
this case, the frequency shifts AJ^ of the antisymmetric 
states are found to be 

Ai = Qf, 

Al = inF + nN)/2-ujB/2, 

Al = inF + nN)/2 + ujB/2, (46) 

where the Bohr frequency is given by 

CUB - V46^ + {nF-nN)^. (47) 

A plot of the frequency shifts A^ as a function of the 
interatomic distance R and for 9 = 7r/2 is shown in 
Fig-IHd). Note that the degeneracy and the level crossing 
of the eigenvalues is removed for 6 [see Sec. lIV A] . 
The eigenstates that correspond to the frequency shifts 
in Eq. (gS]) read 

l<Pa) = 1^2) , 

= e^«sinz?,|V'+) + cos^JV->. 

= -e'^cosda\^i)+smda\i^-), (48) 

where S = \6\e^^ (^ £ {0,7r}), the states Itpf) are defined 
in Eq. ([BS]) . and the angle -da is determined by 

tan2z9a= ^ , < t?, < J . (49) 

I If — ^ 

If the distance between the atoms is small such that 
R < 0.63 X Ao, we have flp < ^n- In this case, we find 
lim^^o \ fa) = |'0a) ^nd lim5_^o A^ = A^, where the eigen- 
states Ii/Iq) and the frequency shifts A^ of the degenerate 
system are defined in Eqs. l|35|) and l|37)). respectively. 

The matrix representation of Ha + -fff2 with respect to 
the symmetric states IV's)} defined in Eq. (1551) 



is found to be 

(uja — ^F S cos 9 \ 
6 cos 9 ujo — ^F —Ssm9 . 
—Ssm9 luq ^ I 

(50) 

Just as in the case of the antisymmetric states, the eigen- 
values of [Ha + Hii]s are written as 

El = h{ua + Al) , 
El = n (c.0 + A^) , 
El = h{uo + A3) , (51) 

and the frequency shifts A* depend only on the inter- 
atomic distance R and the azimuthal angle 9. 

If the atoms are aligned in the x- y- plane [9 — 7r/2), 
the frequency shifts A* of the symmetric states are given 

by 

Al = -^F, 

Al = -{VLf + ^n)I2 + ujbI2, 

Kl = -{nF + nN)/2~UB/2, (52) 

and the corresponding eigenstates are 

\lpI) = -e'«cos79s|V'+) +sini?^|V'7) , 

\^l) = e*«sini?,|V^+)+cosi?,|^7). (53) 

The states ) are defined in Eq. (gO]), (5 = |(5|e'« (^ £ 
{0, tt}), and the angle ds is determined by 

tan2^, = — ^J— , 0<i9, <J. (54) 

5 l]S — ilp ^ 

For small values of the interatomic distance R such that 
^If < ^N, we find limi-^o — IV'I) and lima^o A!, = 
A*, where the eigenstates and the frequency shifts 
A* of the degenerate system are defined in Eqs. ((39)) and 
dH]), respectively. 

Finally, we note that the ground state |4,4) and the 
excited states G {1,2,3}) are eigenstates of 

Ha + Hn- These states together with the symmetric 
and antisymmetric eigenstates of Ha + H^ form the new 
basis of the total state space Tisys- The complete level 
scheme of the non-degenerate system is shown in Fig. [T] 

V. POPULATION OF THE DECOHERENCE 
FREE SUBSPACE 

In this section we describe a method that allows to 
populate the subspace A spanned by the antisymmet- 
ric states. For simplicity, we restrict the analysis to the 
degenerate system [5 = 0) and show how the states 1-00) 
can be populated selectively by means of an external laser 
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FIG. 7: Complete level scheme of the non-degenerate system (5 7^ 0). For the special geometrical setup where the atoms are 
aligned in the x-y-plane (6' — n/2), the analytical expressions for the states \iPa), lifl) and the frequency shifts AJ,, are given 
in Eqs. (|48p . (|53p . (|46p and (|52p . respectively. The frequency shifts AJj (A^) of the antisymmetric (symmetric) states and the 
splitting of the excited states are not to scale. Note that the frequency shifts A^ and A^ depend on the relative position of the 
atoms. 



field. However, a laser field cannot induce direct transi- 
tions between the ground state |4,4) and \tpD as long as 
the electric field at the position of atom 1 is identical to 
the field at the location of atom 2. By contrast, a direct 
driving of the antisymmetric states is possible provided 
that one can realize a field gradient between the positions 
of the two atoms. Since we consider an interatomic spac- 
ing R that is smaller than Ao/2 such that the states in 
A are subradiant, the realization of this field gradient is 
an experimentally challenging task. Several authors pro- 
posed a setup where the atoms are placed symmetrically 
around the node of a standing light field [H, , and this 
method also allows to address the states of our system 
individually. Other methods [H, [s^, El] rest on the as- 
sumption that the atoms are non-identical and cannot be 
applied to our system comprised of two identical atoms. 
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TABLE I: Polarization of the external laser field that couples 
an antisymmetric state IV'a) to an excited state |i, j) £ 
{1,2,3}) for (5 = 0. Note that |?/>^) does not couple to z- 
polarized light, |?/)^) does not couple to y-polarized light and 
I'l/'o) does not couple to x-polarized light. See also Fig. |S1 
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FIG. 8: (Color online) Laser-induced coupling of l^/i^) to the 
excited states |i, j) (i, j £ {1,2,3}) in the case of the degen- 
erate system. States that are not directly coupled to |^^) 
have been omitted (except for the ground state). The laser 
polarization that couples the antisymmetric state to a 
state G {1,2,3}) is indicated next to the respective 

transition. \ipa) is completely decoupled from a j/-polarized 
laser field. 



Here we describe a method that allows to populate 
the states individually and that does not require a 
field gradient between the positions of the two atoms. 
It rests on a finite distance between the atoms and ex- 
ploits the fact that the antisymmetric states may be pop- 
ulated by spontaneous emission from the excited states 
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ihj G {1:2,3}). For a given geometrical setup, we 
choose a coordinate system where the unit vector e^. co- 
incides with the separation vector R. In this case, we 
have 9 — tt/2 and 4> — 0. The ^-direction is determined 
by the external magnetic field and can be chosen in any 
direction perpendicular to R. The polarization vector of 
the laser field propagating in z-direction lies in the x-y- 
plane and can be adjusted as needed, see Eq. (fT7)) . In the 
presence of the laser, the atomic evolution is governed by 
the master equation ^5]) . We find that the coupling of 
the states \ipl) to the excited states \i,j) G {1, 2, 3}) 
depends on the polarization of the laser field (see Table U 
and Fig. [8]). In particular, it is found that \ipl) does not 
couple to 0-polarized hght, \^-^) does not couple to y- 
polarized light and does not couple to x-polarized 
light. At the same time, the states \4>a) are populated by 
spontaneous emission from the excited states. This fact 
together with the polarization dependent coupling of the 
antisymmetric states allows to populate the states 
selectively. In order to populate state for example, 
one has to shine in a y-polarized field. Since the spon- 
taneous decay of IV'a) is slow and since IV'a) is decoupled 
from the laser, population can accumulate in this state. 
On the other hand, the states |'(/'a) and j-i/'a) are depop- 
ulated by the laser coupling to the excited states. This 
situation is shown in Fig. [HUa) for two different values of 
the interatomic distance R. The initial state at i = is 
|4, 4), and for i • 7 = 20 the population of j?/'^) is approx- 
imately 1/4. Since all coherences between ji/;^) and any 
other state are zero, the probability to find the system at 
t — 20/7 in the pure state \tp^) is given by 1/4. 

The exact steady state solution of Eq. ((T8|) is difficult 
to obtain analytically. However, one can determine the 
steady state value of (V'alfi'l'/'a) with the help of Eq. 



hm Clit) /{2Tl) 

t — >-oo 



(55) 



The population of lip^) in steady state is thus hmited 
by the population of the relevant excited states that are 
populated by the y-polarized laser field and that decay 
spontaneously to |V'a)- Furthermore, it is possible to gain 
some insight into the time evolution of (V'alfi'lV'a)- F^i' a 
strong laser field and for a small value of R, reaches 
the steady state on a timescale that is fast as compared 
to l/(2r^). We may thus replace by its steady state 
value in Eq. The solution of this differential equa- 

tion is 



hm cm 

t — >-oo 



2r^„ 



(56) 



and reproduces the exact time evolution of (V'a|p|'0a) ac- 
cording to Fig. UKa) quite well. Moreover, it becomes 
now clear why it takes longer until the population of 
Itp"^) reaches its steady state if the interatomic distance 
R is reduced since the decay rate 2r^ approaches zero as 
R^O. 

So far, we considered only the population of \ipa)i but 
the treatment of and \4>a) is completely analogous. 
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FIG. 9: Time-dependent population of the states \tpa) for 
different polarizations of the driving field. The initial state at 
t = is |4, 4). The parameters are 6 — it/2, (j> = 0, S — and 
A2 = 0. (a) Population of lipa) for Q,y{ri) — Q,y{r2) — 67. 
The states \ipi) and Itpa) are not populated, (b) Population 
of for n^n) = n^ri) = 57. The states and 
are not populated. 



The population of \ip^) by a x-polarized field is shown 
in Fig. [5)^b). The differences between plot (a) and (b) 
arise since the decay rates of \^p'^) and \ipa) are different 
for the same value of R (see Sec. lIVBp . In general, the 
presented method may also be employed to populate the 
antisymmetric states of the non-degenerate system selec- 
tively. In this case, the polarization of the field needed 
to populate a state \lpI} is a function of the detuning S. 

In conclusion, the discussed method allows to popu- 
late the antisymmetric states selectively, provided that 
the interatomic distance is larger than zero. If the in- 
teratomic distance is reduced, a longer interaction time 
with the laser field is required to reach the maximal value 
of (V'al^?IV'a) ~ 1/4. Note that a finite distance between 
the atoms is also required in the case of other schemes 
where the atoms are placed symmetrically around the 
node of a standing light field 0, [1^. While the lat- 
ter method allows, at least in principle, for a complete 
population transfer to the antisymmetric states, its ex- 
perimental realization is difficult for two nearby atoms. 
By contrast, our scheme does not require a field gradient 
between the atoms and is thus easier to implement. It 
has been pointed out that the population transfer to the 
antisymmetric states is limited by the population of the 
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excited states that spontaneously decay to an antisym- 
metric state IV'a)- Although this limit is difficult to over- 
come, an improvement can be achieved if the fluorescence 
intensity is observed while the atom is irradiated by the 
laser. As soon as the system decays into one of the states 
IV'a), the fluorescence signal is interrupted for a time pe- 
riod that is on the order of l/(2r^) (see Sec. UVBt . The 
dark periods in the fluorescence signal reveal thus the 
spontaneous emission events that lead to the population 
of one of the antisymmetric states. 



VI. INDUCING DYNAMICS WITHIN THE 
SUBSPACE A 

In this Section we assume that the system has been 
prepared in the antisymmetric state IV'a), for example 
by one of the methods described in Sec. |Vl The aim 
is to induce a controlled dynamics in the subspace A of 
the antisymmetric states. We suppose that the atoms 
are aligned along the a;-axis, i.e. 9 — tt/2 and (p — 0. 
According to Eq. (|44p . the state li/'a) is then only coupled 
to IV'a)- Apart from a constant, the Hamiltonian Hq 
that governs the unitary time evolution in the space Q 
spanned by {IV'a), IV'a)} can be written as 
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FIG. 10: (Color online) Bloch sphere representation of the 
system dynamics in the subspace Q spanned by the states 
{IV'a), IV'a)}- At t — 0, the system is in the pure state |Va) 
and a static magnetic field is switched on. The Bloch vector 
is rotated around an axis in the a:;-z-plane, and the tilt of this 
axis in a;-direction increases with the magnetic field strength. 
The value of the parameter S is (a) 5 = 3.15 x 7, (b) 5 = 
4.83 X 7 and (c) S = 6.22 x 7, and we chose R = 0.1 x Ao. 



Ho = h 



-5 



nF)/2 -5 

(f7jv - nF)/2 



= fkUBn ■ cr/2 . 



(57) 



where the vector cr = {ax, (Jy, az} consists of the Pauli 
matrices (Ti, and the unit vector n is defined as 



i2S,o,nN-nF)/ujB- 



(58) 



The Bohr frequency lub is the difference between the 
eigenvalues of Hq and is given in Eq. (|i71) of Sec. IIVCI 
Equation (j57p implies that the parameter S which can be 
adjusted by means of the external magnetic field intro- 
duces a coupling between the states IV'a) and IV'a)- If the 
initial state is IV'a), the final state |V'f) reads 



\M))=Uit,OM) 



(59) 



where U = exjp{—iHQt/h) is the time evolution operator. 
The time evolution induced by Hq can be described in 
a simple way in the Bloch sphere picture Q. The Bloch 
vector of the state |V'f(^)) is defined as 



B{t) = (V'F(t)k|VF(<)) . 



(60) 



Initially, this vector points into the positive z-direction. 
The time evolution operator U rotates this vector on the 
Bloch sphere around the axis h by an angle u)Bt- Accord- 
ing to Eq. the axis of rotation lies in the x-z-plane 
and its orientation depends on the parameter 5 which can 
be controlled by means of the magnetic field. In order to 



demonstrate these analytical considerations, we numer- 
ically integrate the master equation ([7]) with the initial 
condition Q{t = 0) = |V'a)(V'al- We define a projector 
onto the space spanned by {IV'a), IV'a)}, 

P=\^l){^l\ + \i>l){^l\- (61) 
The generalized Bloch vector is then defined as 



BN{t)= Tr (TPg{t)P 



(62) 



In contrast to B, Bm is not necessarily a unit vector, 
but its length can be smaller than unity due to sponta- 
neous emission from IV'a) and IV'a) to the ground state. 
Figure [TO] shows the evolution of B^v for different values 
of the parameter 5 which depends on the magnetic field 
strength. Let S — {Sx, Sy, Sz} be a point on the Bloch 
sphere that lies not in the y-z-plane [Sx 7^ 0). If one 
chooses the parameter 5 according to 



2\Sx\ 



\nF-nN\ Sign(5,), 



(63) 



then S lies on the orbit of the rotating Bloch vector 
B if spontaneous emission is negligible. According to 
Eq. (|63p . any point close to the y-z-plane requires large 
values of 5 since \5\ diverges for S'^, — > 0. The dynamics 
that can be induced by a static magnetic field is thus re- 
stricted, particularly because we are only considering the 
regime of the linear Zeeman effect. 

These limitations can be overcome if a radio-frequency 
(RF) field is applied instead of a static magnetic field. If 
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the RF field oscillates along the z-axis, the Hamiltonian 
Ha in Eq. Q has to be replaced by 

Hiit) ^fuJoj^T. ^1+^1- + ^rf(t) , (64) 

i—l fi—l 

where 

2 

Kf(t) - 2M(t) {4^4''- ~ S['is[''}^ (65) 

describes the interaction with the RF field and 

d{t) = do cos(cJrf i + 0rf) ■ (66) 

In this equation, the magnitude of 5o(> 0) depends on 
the amplitude of the RF field, and Wrf and 0rf are the 
frequency and phase of the RF field, respectively. We as- 
sume that the interatomic distance of the atoms is smaller 
than R — 0.63 x Aq. In this case, the dipole-dipole in- 
teraction raises the energy of with respect to 1^^), 
and the frequency difference between these two states is 
fiAT — > 0. Furthermore, we suppose that the de- 
tuning Aif — ujii — (r^Af — Qf) of the RF field with the 
IV'a) ^ iV'a) transition and the parameter Sq are small as 
compared to {rtj^ — Vlp) such that the rotating- wave ap- 
proximation can be employed. In a frame rotating with 
Wrf, the system dynamics in the subspace Q spanned by 
{IV'a); IV'a)} is then governed by the Hamiltonian 

/ Arf/2 -5oexp(i(?!)rf) \ 

Hq = fi\ I 
V -So exp(-i(/)rf) -A,.f/2 / 

^hn,{n,i-(T/2, (67) 

where 

rirf — (— 2(5o cos 0i.f, —2^0 sin ^rf , A^f) /fJrf (68) 

and fly-f = ^A^j -1- 4|(5op. For a resonant RF field (Arf = 
0), the axis riif lies in the x-y-plane of the Bloch sphere, 
and its orientation can be adjusted at will by the phase 
(/)rf of the RF field. Any single-qubit operation can thus 
be realized by a sequence of suitable RF pulses Q. In 
particular, a complete transfer of population from 
to lip^) can be achieved by a resonant RF pulse with a 
duration of vr/firf and an arbitrary phase 0i.f. 

Next we demonstrate that the Hamiltonian Hq in 
Eq. (p7)) describes the system dynamics quite well if the 
atoms are close to each other such that spontaneous emis- 
sion is strongly suppressed. For this, we transform the 
master equation ([7]) with H^J instead of Ha in a frame ro- 
tating with Wrf. The resulting equation is integrated nu- 
merically without making the rotating-wave approxima- 
tion. We suppose that the system is initially in the state 
\4'a)j 3.nd the phase of the resonant RF field has been set 
to (j)ri = T^- Figure [11] shows the time evolution of the 
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FIG. 11: (Color online) Complete population transfer from 
[ipa) to IV'a) by means of a resonant RF field. At t = 0, 
the Bloch vector Bjv points into the positive z direction. At 
t — Tv/Qri, the state of the system is \ipa) and Sjv points into 
the negative z direction. Note that the length of Bat is slightly 
smaller than unity for t > due to the small probability of 
spontaneous emission to the ground state. The parameters 
are R = 0.05 x Aq, 5o = 7, (/"rf = tt and Arf — 0. 



Bloch vector B^- As predicted by Eq. ([57]) . the Bloch 
vector is rotated around the x-axis and at t — n/fl,-{, Bn 
points in the negative z-direction. Due to the small prob- 
ability of spontaneous emission to the ground state, the 
length of Bat is slightly smaller than unity (|i?Ar| ~ 0.95) 
at t — 7r/f2i-f- 

Finally, we briefly discuss how the final state \ipF{t)) 
could be measured. In principle, one can exploit the 
polarization-dependent coupling of the states \ip'^) and 
iV'a) to the excited states (see Sec.|Vl. For example, one 
could ionize the system in a two-step process, where \^p'^) 
(IV'q)) is first resonantly coupled to the excited states 
KjJ) ihj £ {lj2,3}). A second laser then ionizes the 
system, and the ionization rate is a measure for the popu- 
lation of state jf/'a) (iV'a))- Another possibility is to shine 
in a single laser whose frequency is just high enough to 
ionize the system starting from iV'a). Since the energy of 
IV'a) is higher than those of |'0a)j the ionization rate is a 
measure for the population of state \tp^). 



VII. ENTANGLEMENT OF THE COLLECTIVE 
TWO-ATOM STATES 

In Sec. IIV Al we determined the collective two-atom 
states lipl) and that are formed by the coherent part 
of the dipole-dipole interaction. Here we show that these 
states are entangled, i.e. they cannot be written as a sin- 
gle tensor product (g) \%p2) of two single-atom states. 
In order to quantify the degree of entanglement, we cal- 
culate the concurrence [H, of the pure states \^pl} 
and {ipl). The concurrence for a pure state IV'12) of the 
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two-atom state space Tisys = 'Hi ® 7^2 is defined as [i^l 

= V'2[l-Tr(g2)]. (69) 

Here qi = Tr2(^?) denotes the reduced density operator 
of atom 1. The concurrence C of a maximally entangled 
state in Tigys is Cmax = ^73/2, and C is zero for product 
states [131 ■ '^'^d the antisymmetric and symmet- 
ric states \'4>a) and are entangled, but the degree of 
entanglement is not maximal, 

C{m) = C{Wl)) = l<C^... (70) 

Next we compare this result to the corresponding results 
for a pair of interacting two-level systems with ground 
state \g) and excited state |e). In this case, the exch ang e 
interaction gives rise to the entangled states [H, [s^, HI] 

|±) = i=(|e,g)±|g,e)) (71) 

with C(|±)) = 1. It follows that the degree of entangle- 
ment of the states |±) is the same than the degree of en- 
tanglement of the symmetric and antisymmetric states of 
two four-level systems. On the other hand, the states |±) 
are maximally entangled in the state space of two two- 
level systems. This is in contrast to the states j-f/'a) and 
which are not maximally entangled in the state space 
of two four-level atoms. Note that the system of two four- 
level atoms shown in Fig. [1] may be reduced to a pair of 
two-level systems if the atoms are aligned along the z- 
axis. For this particular setup, all cross-coupling terms 
and Fy with i ^ j vanish [see Eqs. and (dH)] 
such that an arbitrary sublevel of the Pi triplet and the 
ground state 5*0 form an effective two-level system. 

In Sec. lVIi we showed that a static magnetic or RF field 
can induce a controlled dynamics between the states ji/'a) 
and I'i/'a)- We find that the degree of entanglement of an 
arbitrary superposition state 

\i^.^^) = a\i^l) + hm (72) 

with |ap + |6p = 1 is given by C(|V'sup)) ~ 1- It follows 
that the degree of entanglement is not influenced by the 
induced dynamics between the states jV'a) and IV'a)- 

Finally, we point out that the antisymmetric states 
IV'a) can be populated selectively, for example by the 



method introduced in Sec. [V] Since the spontaneous de- 
cay of the antisymmetric states is suppressed if the inter- 
atomic distance is small as compared to mean transition 
wavelength Aq, we have shown that the system can be 
prepared in long-lived entangled states. 



VIII. SUMMARY AND DISCUSSION 

We have shown that the state space of two dipole- 
dipole interacting four-level atoms contains a four- 
dimensional decoherence-free subspace (DFS) if the in- 
teratomic distance approaches zero. If the separation of 
the atoms is larger than zero but small as compared to the 
wavelength of the Sq <-*■ Pi transition, the spontaneous 
decay of states within the DFS is suppressed. In addi- 
tion, we have shown that the system dynamics within the 
DFS is closed, i.e., the coherent part of the dipole-dipole 
interaction does not introduce a coupling between states 
of the DFS and states outside of the DFS. 

In the case of degenerate excited states [5 = 0), we 
find that the energy levels depend only on the inter- 
atomic distance R, but not on the angles 9 and (j). This 
result reflects the fact that each atom is modelled by 
complete sets of angular momentum multiplets [s^. We 
identifled two antisymmetric collective states (|'0q) and 
within the DFS that can be employed to represent 
a qubit. The storing times of the qubit state depend on 
the interatomic distance R and can be significantly longer 
than the inverse decay rate of the 5*0 Pi transition. 
Moreover, any single-qubit operation can be realized via 
a sequence of suitable RF pulses. The energy splitting 
between the states |'0a) and IV'a) arises from the coher- 
ent dipole-dipole interaction between the atoms and is 
on the order of IO7 = (10 - 1000) MHz in the relevant 
interatomic distance range. The coupling strength be- 
tween the RF field and the atoms is characterized by the 
parameter (5o which is on the order of (IbBq, where /is 
is the Bohr magneton and Bq is the amplitude of the RF 
field. Since /is is about 3 orders of magnitude larger than 
the nuclear magneton, typical operation times of our sys- 
tem may be significantly shorter than for a nuclear spin 
system. 
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